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^ ' Abstract 

m 

We show that it is possible to determine the locus of Fisher zeroes in the thermodynamic limit for 
f^ I the Ising model on planar ("fat") </!>* random graphs and their dual quadrangulations by matching 

[~^ . up the real part of the high and low temperature branches of the expression for the free energy. The 

f^ ' form of this expression for the free energy also means that series expansion results for the zeroes may 

be obtained with rather less effort than might appear necessary at first sight by simply reverting the 
series expansion of a function g{z) which appears in the solution and taking a logarithm. 

Unlike regular 2D lattices where numerous unphysical critical points exist with non-standard expo- 
2 ■ nents, the Ising model on planar 0* graphs displays only the physical transition at c = exp(— 2/3) — 1/4 

and a mirror transition at c = —1/4 both with KPZ/DDK exponents (a = —1,13 = 1/2,7 = 2). The 
^_^ relation between the </>* locus and that of the dual quadrangulations is akin to that between the 

y , (regular) triangular and honeycomb lattices since there is no self-duality. 

o 
o 
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1 Introduction 

One of the more remarkable results to emerge from the study of various statistical mechanical models 
coupled to two-dimensional quantum gravity is a solution of the Ising model in field [Q, |[ . In discrete form 
the coupling to gravity takes the form of the spin models living on an annealed ensemble of triangulations 
or quadrangulations, or their dual planar graphs. The partition function for the Ising model on a single 
graph G" with n vertices 

is promoted to a partition function which incorporates a sum over some class of graphs {G"} by the 
coupling to gravity 

Z„(/3,/i)= ^ Z,i„g,e(G",/3,/i). (2) 

{G"} 

That introducing such an annealed sum over graphs into the partition function should be a discrete 
version of coupling to gravity becomes clearer when one considers the approach taken to simulating such 
models. In simulations one changes both the geometry (i.e. connectivity) of the lattice and the spins 
on the same timescale, so the spins affect the geometry and vice-versa, mimicking the back-reaction of 
matter and gravitation in the continuum theory. The solution to the Ising model in H] proceeded by first 
forming the grand canonical partition function 

^ = E((r5F)"^"(/3» (3) 

and then noting that this could be expressed as the free energy 

Z = -\ogJv4„Vc^2 exp (-Tr[^{cf>l + q^l)^ccl,icl,2-^{e''4>t + e-''4)]\ (4) 

of a matrix model, where we have written the potential that generates cj)^ graphs. In the above 01^2 
are N x N Hermitian matrices, c — exp(— 2/3) and the TV -^ oo limit is performed in order to pick out 
planar graphs. The graphs of interest are generated as the Feynman diagrams of the "action" in equ.(B|), 
which is constructed so as to weight each edge with the correct Boltzmann weights for nearest neighbour 
interaction Ising spins. Since the edges carry matrix indices the graphs in question are "fat" or ribbon 
graphs. 

The integral of equ.(0) can be evaluated using the results of H to give 
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f" dt , , 1 r dt , ,, 



where g is defined by 



g{z) = ic^z-^ + z 



1 2 6z(cosh/i — 1) 



c2 



(l-3z)2 " ' (l-9z2)2 J- ^^) 

The implicit form of the solution may make it a little difficult to see what is going on, but since 
the singularities of Z w.r.t. g determine the asymptotics of the Z„, the procedure for extracting the 
thermodynamic limit is to look at the solutions to g'{z) = 0. These can be explicitly determined when 
h = Q a.s ZL ~ —1/3, z_f/i.2 = (1/3) [1 T -j=\, ^-^3.4 ~ (1/3) [1 ± 7=^] ^^^ then substituted into g 



(7) 
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where z^ is the low-temperature branch and the various Zfj high-temperature branches. Inserting the 
annronri ate n(r\ into the exnression for Z„ 



appropriate g{c) into the expression for Z„ 

-Acg{c) 



n-" 



_(1-C2)2 

then gives the asymptotics of the Z„ and the thermodynamic behaviour of the free energy per site F 

-4:cg{c 



(8) 



i^ = - log 



(l-c2) 



,2^2 



(9) 



since if we are given a canonical partition function Z„ the associated free energy per site F will be given 
in the thermodynamic limit by 

F^ lim -logZ„. (10) 

n — ^oo 77, 

The third order phase transition with the so-called KPZ/DDK [Q exponents, a = —1, /? = 1/2, 7 = 2 
occurs when ghic) = gHiic) which gives a critical coupling c = 1/4. It is possible to carry out a 
perturbative expansion in h around the h = solutions above to obtain the magnetic critical exponents 
directly from the discretized formulation [^ and it is likewise possible to confirm universality by solving 
the model on (f)^ graphs. The KPZ/DDK exponents were verified in a continuum formalism in y using 
conformal field theory techniques. 

Given the solution of ||l|, g] it is tempting to use it as a test case to investigate various statistical 
mechanical ideas and methods, in much the same manner as the Onsager solution has served as a paradigm 
over many years. One such effort was presented in Iq], where the behaviour of the partition function zeroes 
for the Ising model coupled to two-dimensional quantum gravity was investigated by series expansion and 
numerical means. The study of partition function zeroes for statistical mechanical models was initiated 
by Lee and Yang for complex external fields p, m and later extended by Fisher and others to complex 
temperatures B. It offers an alternative viewpoint of the approach to the thermodynamic limit and 
means of extracting critical exponents. A study of partition function zeroes for the Ising model coupled 
to two-dimensional gravity addresses several interesting questions. It is not clear a priori that loci of 
partition function zeroes will continue to lie on simple curves in the c = exp(— 2/3) or y = exp(— 2/i) 
planes when a sum over some class of graphs, in this case planar graphs, is folded into the partition 
function. Although this is generically the case for the Onsager and related solutions on regular two- 
dimensional lattices fi Iffl, there are exceptions such as the "bathroom-tile" lattice. Other sorts of 
behaviour are possible too. For instance, introducing geometric disorder in the form of Penrose tilings 
gave an complicated extended structure of temperature zeroes away from the physical critical point, but 
still gives rise to Onsager exponents |1^. Fractal lattices on the other hand display an intricate fractal 
pattern of zeroes Jl^ . 

The work in (5| suggested strongly that the temperature zeroes did lie on curves and that the field 
zeroes still lay on the unit circle in the complex y = exp(—2h) plane, as in the regular lattice Onsager 
solution. Similarly, a comprehensive analytical study of the Lee- Yang zeroes for the Ising model on planar 
(j)'^ graphs was carried out in [|l3[ , where it was found implicitly that the Lee- Yang circle theorem still 
held, since the complex field singularities were shown to lie at purely imaginary field values. In this paper 
we concentrate on the temperature (Fisher) zeroes for the Ising model on planar cj)'^ graphs and their dual 
quadrangulations, showing how to derive the locus of zeroes analytically using the idea that the locus 
should be thought of as Stokes lines 1^, |l^, |lj, 14, |5|. We compare the results with the various sorts of 



behaviour observed in ||^, [l^l for regular two-dimensional lattices and also note that the form of the Ising 
solution means that the various zeroes on finite planar (j)'^ graphs can be extracted without evaluating a 
series expansion for the full expression for Z. In the sequel we first briefly discuss the general background 
to Lee- Yang and Fisher zeroes and the analytic determination of the loci of zeroes. We then show how 
series expansion results for finite graphs, such as those in |5), can be recovered and extended economically 
before we move on to discuss obtaining the loci of Fisher zeroes for the Ising model on cj)'^ graphs and 
their duals analytically. 



2 Lee- Yang and Fisher Zeroes 

The starting point of Lee and Yang's work [g, g was the consideration of how the non-analyticity char- 
acteristic of a phase transition appeared from the partition function on finite lattices or graphs, which 
was a polynomial 

Z = Y,D,r,nC^y^ (11) 

for a lattice with m edges and n vertices, again with c — exp(— 2/3), y — exp(— 2/i). They showed that 
the behaviour of the zeroes of this polynomial in the complex z plane, in particular the limiting locus as 
m,n ^ cx), determined the phase structure. Similarly, the behaviour of the zeroes in the complex c plane 
determines the nature of temperature driven transitions |Q. In the latter case, in zero external field for 
simplicity, the free energy on some lattice or graph G„ with n nodes and m edges can be written 



F(G„,/3)^-lnn(c-Cfe(/3)), (12) 

which in the thermodynamic limit becomes 



fe=i 



F{Goo.P) - - / dcp{c)\n{c-c{L)), (13) 

where L is some set of curves, or possibly regions, in the complex c plane on which the zeroes have 
support and p{c) is the density of the zeroes there. The singular behaviour of p(c) as c approaches the 
physical transition point cpT is related to the specific heat exponent a by 

p{c)^{c-cpTf-^. (14) 

The general question of how to extract the locus of zeroes analytically has been considered by various 
authors. It was observed in ||l^ that such loci could be thought of as Stokes lines separating different 
regions of asymptotic behaviour of the partition function in the complex temperature or field planes. 
Across these lines the real part of the free energy should be continuous and the discontinuity in the 
imaginary part should give the density of zeroes. Shrock and collaborators ^, |l^, [ij] have obtained 
many interesting and explicit results on the Fisher and Lee- Yang loci for the Ising and other models 
on regular lattices by matching free energies in this manner. They also observed that the condition 
3?^~^ — 0, where ^ was the correlation length, gave equivalent loci jlj]. Both this condition and the 
matching of free energies are consistent with the idea that the loci of zeroes coincide with a change of 
dominant behaviour in the asymptotics. 

More recently, the case of models with first-order transitions has been investigated by Biskup et al. 
llq] who showed rigorously n that the partition function of a d-dimensional statistical mechanical model 
defined in a periodic volume V = L'^ which depends on some complex parameter such as c or j/ can be 
written in terms of complex functions Fi describing k different phases as 

Z = ^g,e-'='^'^ + 0(e-^/^°e-^^^), (15) 

where qi is the degeneracy of phase Z, /3 is the inverse temperature and io is of the order of the correlation 
length. The various Fi are the metastable free energies per unit volume of the phases, with 5ftF; — F 
characterising the free energy when phase I is stable. The zeroes of the partition function are then 
determined to lie within 0{e~^^^°) of the solutions of the equations 

^F^f = 5RF^^{ < ^F',f/, Vfc ^ I, m, 

PVi^Fi^L - ^F^,l) = TT mod 27r. (16) 



^ Under suitable technical conditions. 



The equations ( pij| ) are thus in perfect agreement with the idea that the loci of zeroes should be Stokes 
lines, since the zeroes of Z asymptotically lie on the complex phase coexistence curves 'SiFij^ = ^F„ii^ 
in the complex parameter plane. 

The specific Biskup et al. results apply to models with first order transitions - the canonical example 
being the field-driven transition for the Ising model, and we are interested here in a model with a third 
order transition, so it might initially seem that these results were inapplicable. We are saved by the 
fact that when considered in the complex temperature plane the transition is continuous only at the 
physical point itself (and possibly some other finite set of points). This is nicely illustrated by looking 
at expressions for the magnetization for the Ising model on the square lattice, on fat (planar) ^^ graphs 
and on thin (generic) (j)^ graphs: 



\i/ 



(1 + u)1/4(i_6u + m2)1 

^ = ^ — (i_^)i/2 — — (^q'^^'-e)' 



3(1 -16u) 



1/2 



M = "^% ^":' (fat 04), (17) 

n _ v^i/2 

" - (1-20(1 + 0./. <"-*'). 

where u = c^ — exp(— 4/3). It is clear from these expressions, which apply through the complex extension 
of the low-temperature phase with M zero outside, that although M will vanish continuously at the 
physical critical points, u = 3 — 2^2; u = 1/16 (i.e. c — 1/4); c — 1/3 respectively ^, it will generally be 
non-zero at the phase boundary approaching from within the low-temperature region, whereas it will be 
zero approaching from outside, which is characteristic of a first-order transition. 

In summary, both general considerations about the change of asymptotic behaviour of expansions of 
the partition function in different regions of the complex temperature or field planes |l3, W\ Bl nO| and 
rigorous results |l£| lead to equ.(|l6[) as a means of determining the loci of zeroes. 



3 Series Expansions With (two thirds) Less Pain 

To get a series expansion for Z one must in principle go back to equ.dq) and invert (or more correctly, 
revert) the expression for g{z) in equ.(o) expanded as a series in z to get an expansion z{g) in powers 
of g. This is then inserted in equ.(p^ in order to obtain the desired series from which the zeroes may 
be extracted. However, if one considers the various terms in cqu.(P) independently some interesting 
observations can immediately be made. Taking each of the terms in equ.(H) separately, 

-^1 : n log I - 
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Zi: — -git), (18) 

9 Jo t 

the series expansion of each term k — 1,2,3 can be written as Zk — J2n'^k^n{c'^)g^ where A„(c^) is 
identical for all the Z^ ■ In addition, normalizing a„ = 1 for the i log I - ) term (since any n dependence 
in this can be put in A„), one finds 

a? = 1, 

77 

"2 - 7 ^, (19) 

2 (n + 2)' "■ ' 

^There are further points where the magnetization vanishes continuously: at the anti-ferromagnetic point m = 3 + 2^2 
and the unphysical point u = —1 on the square lattice; and the unphysical point c = —1/4 on the planar <f>^ graphs, but 
these are discrete and finite in number. 



"^ " (n+l)- 

Although the A„(c^) which determine the partition function zeroes for a given power of g are the same 
for each Zk, this is obscured by the different a]J in the sum Zi + Z2 + Z3. 

Why should this structure be present? The solution given in equ.(||) comes from integrating an 
expression of the form 



Z= I dx{\ 
Jo 



x)\og{f{x)) + ... (20) 



which is common in form to all matrix models. The particular details of a given model are encoded in 
the f{x) which for the Ising model with h = satisfies 



g^- — L. .of.. -^^ +3^M— • (21 



V c ; V(i-¥)' 



V c 



The expression in equ.(ra) emerges on defining z = -|^ and integrating by parts. 

If we expand the log(/(a;)) in the integrand of cqu.(pG), then we get / dx{l~x)P{gx^ c), where P{gx, c) 
is a power-series in gx. Looking at the structure of P{gx, c) it is clear that the integration over x only 
affects numerical factors, but not the c^-polynomials which determine the zeroes. If we now carry out the 
partial integration on the P(gx,c), we obtain three terms: P(l,c)/2 from the boundary corresponding 
to Zi (and irrelevant additional terms); — J dxxP'{gx,c) corresponding to Z2; and ^ ^ dxx^P' {gx,c) 
corresponding to Z3. Carrying out the differentiation of the power series P'{gx,c), followed by the 
integration above recovers the observed values of the ajj . 

The upshot of all of this is that for the purposes of extracting partition function zeroes it is sufficient 
to simply consider the expansion of log(z{g)/g) in powers of g, by reverting the series for g{z), 

5 = (1 - c^)z + 6z^ + 3(c2 + 9)z^ + ... , (22) 

to get 

5(5) = 5 - 652 + 3(c2 + 5)(c2 + 3).g^ + . . . (23) 

and then taking \og{z{g)/g), where we have rescaled z ^ (c^ — l)z, g ^ {c^ — l)^g for algebraic con- 
venience. The polynomial in c^ in front of the appropriate power of g will then yield the desired Fisher 
zeroes. Various efficient algorithms exist for the reversion of series (i.e. getting from equ.(p2|) to equ.(|2^) 
and we have used both the built in algorithms in Maple and Mathematica and one of the earliest numerical 
algorithms, Newton iteration, to revert the series for g{z) |lq], all with identical results. 
For the Newton iteration when /i = we take our starting function to be 

/ = (c2 - l).g(3(c2 - 1)S - 1)2 - 5 + (c^i - 3c2(c2 - I)2i3)(3(c2 _ 1)2 _ 1)2, (24) 

This is then iterated with the standard Newton formula |16{ 

Zk+i = Zk- .- . (25) 



with the starting condition zq = g (see equ.(|23|) above). Since Newton iteration displays quadratic 
convergence the iteration index k is related to the order of the expansion i for z m g hy i = 2''. We 
thus get order 2*^ — 1 for z{g)/g in k iterations. This is both an advantage and a disadvantage since, 
although long series are generated quite rapidly, they are doubling in length at each iteration which can 
rapidly exhaust the available memory. With the built in functions on the other hand, one can proceed 
incrementally in the order. 

To verify that \og{z{g)/g) really is sufficient to determine the zeroes correctly we can compare the 
results for the zeroes in c coming from the polynomial coefficient at a given order in the expansion of 
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log{z{g)/g), for instance g 
84768120 



7 



5255623440c^^ + 547079930820c^ 



17774272305360c2^ + 2565884409300000^" + 1908495144456480c 



18 



7988644803377340c' 



138652618561302240 
c 

7 



30012882991193160c' 



28216084061998800c^° + 16541610886750140c^ + 60022315957168800^ 



1335148577661600c* + 173901 100089600c^ 



95938227092700 

7 ' 



(26) 



with the results from the full expression for Z from m. n There is complete agreement between the 
numerical values of the zeroes obtained with either method as shown in Table 1, 



17.26983082/ 


9.620359803/ 


4.237939134/ 


3.307585457/ 


2.153341531/ 


1.952696297/ 


±.2259213695+ .3562012608/ 


±.1989588142+ .6083974700/ 


±.1697220421 + .9027390050/ 




±.1285487771 + 1.322771774/ 





Table 1: The zeroes from an expansion of both \og{z{g)/g) and Z to order g^^ are identical. 
The complex conjugates of all the values shown are also zeroes. 

It is easy to obtain an expansion of log{z{g)/g) in g up to quite high order with relatively modest 
computing facilities. The coefficient of g'^^ from such an expansion is given in Appendix A and we use 
this for comparison with the analytical expressions for the loci of Fisher zeroes in the next section. 

The observations above regarding the sufficiency of \og{z{g)/g) for determining the partition function 
zeroes also apply to both Lee- Yang zeroes and the Fisher zeroes in non-zero field, since the general 
structure of the expression for Z — J^ dx{l — x) \og{f{x)) + ... is unchanged when the field is turned on 
- it is the defining equation for f{x) which is altered. A nice confirmation of this can be obtained by 
using the series expansion of log{z{g)/g) to obtain the field zeroes in the variable y = exp(— 2/i) which 
are plotted in Fig. 1 for an expansion up to 0(31) with c = 1/4. These clearly lie on the unit circle, as 
they do for the full partition function, and are also evenly distributed. Similarly, Fisher zeroes for the 
partition function in field can also be investigated by expanding \og{z(g)/g), using the full expression in 
equ.(|6|) with h ^ 0. The flow observed reproduces that seen in Q where the complete partition function 
was considered (at much lower order). 



4 The Locus of Zeroes on 0^ Graphs 

The determination of the locus of Fisher zeroes in the thermodynamic limit turns out to be rather 
straightforward, as we now describe. Since we wish to match ^F between the various solution branches 
to obtain the loci of Fisher zeroes and from equ.(10h F ~ log(5(c)), the equation which determines the 
loci of zeroes in the thermodynamic limit is 



or more concisely 



log|gL(c)| =log|5ff,(c)|. 



\9Lic)\ = \gHM\ 



(27) 



(28) 



where the various g are given in equ.(y) and i — 1,2,3,4 where appropriate depending on the region of 
the complex c plane. The explicit expressions arising from substituting a complex value of c into equ.(|2q) 

^We would like to thank the authors of [BJ for providing us with their original data. 



are not very illuminating and we do not reproduce them here, but they allow the resulting curves to be 
plotted with Maple or Mathematica. 

For comparison it is useful to refer back to the locus of Fisher zeroes for the Ising model on a regular 
square lattice, which is the limagon in the complex u = c^ = exp(— 4/3) plane shown in Fig. 2. The use of 
u has the advantage of subsuming the c —^ —c symmetry that is present in the solution and is perhaps the 
most natural choice of variable. The (complex extended) ferromagnetic phase lies inside the inner loop, 
the paramagnetic phase between the loops and the antiferromagnetic phase in the exterior. The physical 
ferromagnetic and antiferromagnetic transition points lie on the positive real axis at m = 3 =F v2 and a 
multiple point with non-standard exponents is present at u = — 1, as already noted in the introduction 
when discussing the magnetization. The limaQon maps onto a pair of overlapping circles in the complex 
c plane, which is probably a more familiar presentation. 

In contrast only two phases are present in the diagram for the Ising model on planar cj)'^ graphs in 
the u plane, since there is no antiferromagnetic phase in this case (the graphs are not loosely packed 

- both odd and even loops can be present). The locus of Fisher zeroes in the u plane for the Ising 
model on planar 0"* graphs is shown in Fig. 3. The interior of the loop is the ferromagnetic phase 
and the exterior the paramagnetic, with the physical transition lying at w = 1/16. The cusp point at 

— j{-^+5\/6) = —6.124 362 . . . does noi represent an unphysical phase transition point, unlike the multiple 
point on the limagon. as can be confirmed by looking at the discriminant A{g') oi eg' (z) {1 — z)^ / (1 — c^)"^ . 
This will show up any non-generic points where multiple roots exist giving phase transition points as 
opposed to the generic first-order boundaries Q. The discriminant is proportional to 



{1-u) 



M9) oc " A 'Z:' (29) 



so the only place where non-generic behaviour appears apart from the trivial point, u = 0, is at the 
transition point u = 1/16. At this point the KPZ/DDK exponents a — —1, /3 = 1/2, 7 = 2 are 
manifested rather than the flat lattice Onsager exponents. In addition to the plotted curve a cut runs 
down the negative u axis from the cusp to w = — oo. In this respect the locus of zeroes is more similar to 
that of the Ising model on a regular triangular lattice rather than the square lattice |^ . Obtaining zeroes 
for finite size graphs from the series expansion one also finds zeroes lying along the negative u axis. 

Back in the complex c = exp(— 2/3) plane one has the locus shown in Fig. 4 with a physical transition 
point at c = 1/4 and a mirror image at c = —1/4, both with KPZ/DDK exponents. We have again omitted 
the cuts running up and down the imaginary axis from the cusp points at ±|(|-f \/6) — ±7 2.474 744 872... 
to ±Joo for clarity. In the c variable it is clear that these cuts, and the associated zeroes, appear because 
of the y/c terms in the paramagnetic solutions gjji in equ.(p[). Due to the cut on the imaginary axis the 
entire left hand region of the c plane exterior to the locus plotted in Fig. 4 forms an unphysical phase 
(labelled "O" in the parlance of pj ) , whereas the exterior right hand region forms the (complex extended) 
paramagnetic phase. The O phase vanishes in the u plane since the sections of the cut on the imaginary 
axis in the c plane are folded onto the negative u axis. Nonetheless the presence of the cut means that 
it is not possible to make a circuit of the origin in the u plane in the paramagnetic region. 

We have also plotted the numerically determined roots obtained from the coefficient of g'^^ in the 
expansion of \og{z{g)/g) in Figs. 3,4. It is clear that the agreement with the analytically derived locus 
is very good for the roots closest to the critical point (s). The finite-size effects, which we have not 
investigated in detail, would be interesting to explore further since it can be seen that in the u plane the 
roots initially lie outside the analytic locus before moving inside it as the cusp is approached. 

The Fisher zero locus on 0'* graphs is thus, if anything, simpler than the square lattice locus - the 
antiferromagnetic phase is absent and there are no unphysical transition points such as the u = —1 point 
on the square lattice. The topology is closer to the locus of zeroes observed for regular triangular lattices 
which are also not self-dual, though the cut(s) do not extend into the ferromagnetic region as they do on 
the triangular lattice. 



*The factor of c/(l — c'^)^ has been included to match with the factors in the free energy in equ.(y) and the (1 — 2)^ 
cancels the (irrelevant) denominator of g'(z). 



5 Dual Quadrangulations 

Although the Ising model on planar 0^ graphs does not display an antiferromagnetic transition it does on 
their dual random quadrangulations |1^ . Again it is useful for orientational purposes to look at the locus 
of zeroes for the square lattice Ising model which corresponds to the case of regular quadrangulations. 
This gives the overlapping circles shown in Fig. 5. 

The locus of zeroes for the random quadrangulations can formally be obtained by mapping the locus 
of Fig. 4 in the c plane to the c* plane. 

which may also be written as c* ~ tanh(/3). The identification of the phases is, however, different for 
the dualized 0'* diagram and the quadrangulations. In Fig. 6 we can see that this gives a ferromagnetic 
phase around the origin, a crescent shaped paramagnetic phase and an exterior region forming the 
antiferromagnetic phase. The boundaries of the crescent are inverse to each other in the unit circle, a 
section of which also forms the image of the cut along the imaginary axis in Fig. 4 and joins the horns 
of the crescent to c = -1 to form the FM/AFM boundary. These horns lie at -0.719 273 115 4... ± 
70.694 727418 2... and just as for the cusp point(s) in the u and c planes there is no sign of non-generic 
behaviour there. The numerically determined roots have again been plotted for comparison. 

One can also obtain the locus for the quadrangulations directly by using the dualized expressions for 
g, replacing c ^ (1 — c)/(l + c), 

1 5c2-22c + 5 
9L ■' 

9H1.2 = ~i— — n"^^-,^2'^"'^ — -^ (31) 

5^3.4 

and plotting |.gL(c)| — \gH{c)\- The results are identical whichever method is used. The ghic) is now 
appropriate for the paramagnetic phase in Fig. 6, and the gnicj^s cover the ferromagnetic and antiferro- 
magnetic phases. 

The dual quadrangulations display a ferromagnetic transition at c = 3/5 (the image of c = 1/4) and 
an antiferromagnetic transition at c = 5/3 (the image of c = —1/4), both with KPZ/DDK exponents. 
Looking at the duahzed discriminant (c.f equ.(|29|)), 

~ (3c -5)^(5c- 3)^(1 + c)i4(l-c)i^ 

Aoc ^jg , (32) 

we can see that this also vanishes at c — —1 as well as the transition points at 3/5, 5/3. It would be an 
interesting exercise to determine the associated exponents in full from the matrix model formulation or 
series expansions, but we do not pursue this further here. 

6 Discussion and Desiderata 

We have seen that it is possible to obtain the locus of Fisher zeroes for the Ising model on both planar 0* 
graphs and their dual random quadrangulations by simply matching the moduli of the critical couplings 
g{c) or g(c) for the various phases. This gives the loci in Figs. 3/4 and 6 respectively which form an 
interesting contrast to their square lattice counterparts in Figs. 2 and 5. In particular, the lack of self- 
duality on the (/)^ graphs makes for less structure than on square lattices since there is no antiferromagnetic 
transition and the general picture is more similar to the loci seen for regular triangular lattices and the 
dual honeycomb lattice g. 

Although the considerations of ||l^ were for first-order transitions the fact that the phase boundaries 
of the complex extended phases generically display first-order behaviour means that the idea of matching 
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up the real part of different branches of the free energy also applies to the models here. The more 
general approach of ||l2|, Q and the idea of regarding the loci of zeroes as Stokes lines or complex phase 
boundaries also suggests that this will determine the loci correctly. Comparison with series expansions 
results, both from B and the much longer self-generated series here, found good agreement with the 
analytically determined loci. We saw that generating such series expansions to obtain either the Fisher 
or Lee- Yang zeroes for the Ising model on planar graphs reduced to the reversion of a series for g{z) and 
taking a logarithm, both simple operations, and relatively long series could be produced without undue 
effort. 

It is perhaps worth noting that various properties of Fisher zero loci discussed in pO| which were 
derived from general considerations still apply to the models considered here. We have: 

• The loci of points where the free energy is non-analytic is symmetric about the 5Rc axis ( "Theorem 
1" of flQ]) - applies to both 0^ graphs and quadrangulations. 

• If the graph has even co-ordination number the locus of zeroes is invariant under c — > — c ( "Theorem 
2" of Q) ~ applies to the (f)^ graphs. 



• If the graph is bipartite the locus of zeroes is invariant under c -^ 1/c ( "Theorem, 3" of |10 ) - 
applies to the random quadrangulations. 

We have looked only cursorily at Lee- Yang zeroes in this paper, for which a circle theorem in the 
complex y = exp(— 2/i) plane still appears to be valid on both planar [^ and thin graphs Q. We noted 
that the results for the series expansions still hold when ft. 7^ 0, and a series in y to obtain the Lee- Yang 
zeroes could still be generated by reversion oi g{z). Since the form of Z is identical for (f)^ planar graphs, 
an investigation of Fisher zeroes on planar (jfi graphs and their dual triangulations along the lines of the 
work here would thus be perfectly feasible. This would make for an interesting comparison with other 
regular lattices and the cases discussed here since, as we have seen, the structure of the loci of Fisher 
zeroes is highly non-universal even when the phase transition in the Ising model is universal. 

Finally, we note that similar methods to those employed here may be used to obtain the locus of 
Fisher zeroes for the (mean field) Ising and Potts model on thin random graphs |19| . 
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8 Appendix A 

The polynomial in u which appears as the coefBcient of g^^ in the expansion of \og{z{g)/g). This was 
used to determine the zeroes plotted in Fig. 3 for comparison with the analytically determined limiting 
locus. 

1560 ti''^ 
19586720 u^^ 



- 37139835653499716006137093004620701844717S3759SS605 

- 14810541168502558092136894448394S303066640247852769jou,^u ^ 

- 219796631965060385324004998785213226455140375358517705237040 u''^ 

- 164831966491059030819403027184073178546645417896670585118583360 u"^^ 

- 73190682361368063761602209974276459964883309418605799340499616080 u''"' 

- 21193244550424189508255837901963813871588047644811792548846191056480 u"^^ 

- 4271156254827515627988674042672332480042287504111958952442525040283760 u^^ 

- 627862 14739052064 1370446437291 1S33507904S0810579 19073 114702243 71 19076480 u"^^ 

- 69745140248196358151864786123993807529941628810020280214779959586962178000 ^l^" 

- 6018464875673457244143786208765208352061116125920486582535296176941628552480 ^l^^ 

- 412471197424814031713421607686052888384939032979910093299169932719217502605680 u^ 

- 22861804785909511304850003848078876423860505468001445816621100351727090582759360 

- 10403880010095256270040385 117928530 1449941 56334662427017510227887386142237818692C 

- 3937229409789412050834300775784434139478345858309842059514S520582863335631694924e 

- 12526699796982776 171 12 1043828330714902967574185029094471295802084671586332292926E 

- 3382423594 17628093 16986439403557899 11983064362 180679589508369104 1491223423768922S 

- 78150699337533941962786903578124700173639477718261805909729851999256785693367615e 

- 155626583739012094142127606161466278147378880000322572229656722886581206971869687 

- 268811911671136955690918042236856880858371066215003407241731389769889134493580751 
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^21 
17 



14 



182 12688770 15292037233 10678 198571 583992864167528 11 996600 134670652 194S13938S 136 78 1691370293 11984688 169520 u^ 
805308804562459495708503916394987237038843253612181001165500116701 1932667983907251509818773391710450560 u^^ 
3252470645082229356008650751768353156846940016343565742435144703608857453346650919121596251905872636880 u^^ 
1199037878489113778438716271817391517512781176088717654560134880626972165041017913820459566740270632480 
403173861504807914325983089651615786813837710424803782422307586554827984880680408117617036296222317680 • 
12354573793 1826333990589263018372737196 158771 1367052 1S30889040933715736389768325488628640605162 1360320 ' 
34469806039566521127126491885908195430460993244436123158311930137387052183055924535551928957561135760 
8747581221336382206763336472263174629514640810049804874785909198519454278642989719373646252543135840 i 
20 1696690407326145679 19067832888779716940 183877445442853869663766449368 11 74326269 15415 10272304280880 i 
4220401 10288334497759627356286466867221 28760259008 14914685925405703837716715472938483132827955 12320 u 
80036988814102662395545702671530110662277119050361283201390143378307264063442392366336074807002080 
13737160482867657888903979612359904387039625719432396863706607709167730140315331444552427440036800 
2130581489746179897063191463661629432483663971226284474690803480226230394996040706276763337037600 ■> 
298084582170998682522557843634130692118118248948169065379102074737434795961756416684382726001600 u 
3754496621550752332222769467284 10 1004372283584497434990642 1354497206822261 6198698856102 12620000 -u?-' 
4247074844108597663669728712049132720278802778418177142466258256291478726641887126832109121600 u^ 
430155384350934174771440197370772245548360430019275647467244387233120149574187173959697949600 u^ 
38846640456271627441559966602880666242219998526979599170426784619478931564006336448368211200 i^ 
3109241968691735189036674466756079403286981061864652920257364770194926524810059222607549600 
218519863192071952043562268682613228805967293880235742328799873278069522182085472663592000 
1328346841532081 16241 19223651866 16482653802801 696412 19315791430662801243870445SSSS3314400 • 
680683119521593923666628073413708660282239032204436877326746488850479986486890462153600 
28058580972988347971863829784009239496026729839677340527392289791671643157102437191200 u' 
845523959563531370690706608441021047213998439269969130238386790640994054828513707200 -u 
1140887404597421591776564792039966147370831359666672145202708266673791740119521072800/79 
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Figure 1: Lee- Yang zeroes in the complex y — exp(— 2ft,) plane for the Ising model on planar cf)^ graphs 
calculated from an expansion of z{g)/g to order y^^ . They clearly lie on the unit circle. 
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Figure 2: The locus of Fisher zeroes in the complex u = exp(— 4/3) plane for the Ising model on the 
square lattice. The ferromagnetic phase lies inside the inner loop, the paramagnetic phase between the 
loops and the antiferromagnetic phase in the exterior. 
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Figure 3: The locus of Fisher zeroes in the complex u = exp(— 4/3) plane for the Ising model on planar 
(j)'^ graphs. Some of the zeroes calculated from the series expansion to order 5^*^ in Appendix A are also 
plotted for comparison, though we have not plotted the zeroes lying at large negative u values on the 
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Figure 4: The locus of Fisher zeroes in the complex c = exp(— 2/3) plane for the Ising model on planar 
(j)'^ graphs. The interior of the curve is the ferromagnetic FM region and the exterior the paramagnetic 
PM and unphysical O phases, separated by the cuts on the imaginary axis which we have not plotted. 



16 



AFM 




-2-1 



Figure 5: The locus of Fisher zeroes in the complex c plane for the Ising model on the square lattice. 
The O region is not connected to any physical values of c. 
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Figure 6: The locus of Fisher zeroes in the complex c plane for the Ising model on random quadrangula- 
tions. This is formally identical to the mapping of the locus in Fig. 5 into the c* = (1 — c)/(l+c) = tanh(/3) 
plane but the identification of the phases is different. For quadrangulations the interior of the crescent 
is the paramagnetic (PM) phase and the two horns of the crescent are joined to c = — 1 by arcs of the 
unit circle which are the images of the imaginary axis cuts in Fig. 4. These form the boundary between 
the inner FM and exterior AFM phases. 
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